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Abstract 


The SPEI (Standardized Precipitation-Evapotranspiration Index) is a variant 
of the widely used SPI (Standardized Precipitation Index). It has greater 
potential as a drought index since it is sensitive to climate change, given that 
it uses a more realistic measure of the water availability, namely, the climate 
balance. This measure is represented by the difference between monthly 
precipitation and potential evapotranspiration. This study describes in detail 
the calculation of the SPEI and its application to the monthly historical 
precipitation and average temperature records available from the Zacatecas 
weather station, located in the capital of the state of Zacatecas, Mexico. The 
data covers 86 years over the period 1930 to 2015. For historical droughts, 
the SPEI results indicate approximate values of 32.2, 11.5, 5.0 and 1.3% for 
light, moderate, severe and extreme droughts, respectively. For this 
particular location, three climatic change scenarios were proposed: the first 
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one adopts a progressive and linear reduction of 20% in annual precipitation; 
the second one adopts a progressive and linear increase of 4 *C in the 
average temperature, and the third superimposes the two effects mentioned, 
which is therefore the most critical. For the three proposed scenarios, the 
results suggest that light and extreme droughts slightly increase in number 
or percentage, and that moderate and severe droughts decrease also only 
slightly. Results are shown with three types of graphs: the SPEI evolution for 
12 months, a historical graph, and third clearly show the times of occurrence 
of drought periods as well as specific differences in the SPEI value. Based on 
these analyses, a systematic application of the SPEI for the characterization 
of historical and future droughts is recommended. Its application for future 
droughts would only be indicative. 


Keywords: SPEI index, Log-Logistic distribution, probability weighted 
moments, potential evapotranspiration, statistical homogeneity, climate 
change. 


Resumen 


El SPEI (Standardized Precipitation-Evapotranspiration Index) es una 
variante del ampliamente difundido SPI (Standardized Precipitation Index), 
que tiene un mayor potencial como índice de sequías, al ser sensitivo al 
cambio climático, debido a que usa una medida más real de la disponibilidad 
de agua: el balance climático. Esta medida se evalúa con la diferencia entre 
precipitación y evapotranspiración potencial mensuales. En este estudio se 
describe con detalle el procedimiento operativo del SPEI y se aplica a los 
registros históricos mensuales disponibles de precipitación y temperatura 
media en la estación climatológica Zacatecas, ubicada en la capital del 
estado de Zacatecas, México, los cuales abarcan 86 años (1930-2015). Para 
las sequías históricas, los resultados del SPEI definen para las sequías leves, 
moderadas, severas y extremas los siguientes valores aproximados: 32.2, 
11.5, 5.0 y 1.3%. Para esta ubicación del país se consideró conveniente 
formular tres escenarios del cambio climático: el primero acepta una 
reducción progresiva y lineal de 20% en la precipitación anual; el segundo 
adopta un aumento progresivo y lineal de 4 *C en la temperatura media 
anual, y el tercero sobrepone los dos efectos citados y, por lo tanto, es el 
más crítico. En los tres escenarios planteados se encuentra que las sequías 
leves y extremas aumentan de forma somera en número o porcentaje, y que 
las sequías moderadas y severas disminuyen también de manera mesurada. 
Se exponen los gráficos de evolución del SPEI de duración de 12 meses, 
histórico y del tercer escenario, en los cuales se visualizan con claridad las 
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épocas de ocurrencia de los periodos de sequía, así como de las diferencias 
puntuales en el valor del SPEI. Con base en los resultados de estos análisis, 
se recomienda la aplicación sistemática del SPEI en la caracterización de las 
sequías históricas y futuras. En estas últimas, los resultados tienen sólo un 
carácter orientativo. 


Palabras clave: índice SPEI, distribución Log-logística, momentos de 
probabilidad ponderada, evapotranspiración potencial, homogeneidad 
estadística, cambio climático. 
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Introduction and objective 


Droughts are one of the main natural causes of damage to agriculture, as 
well as to all economic activities and to the environment. Droughts occur 
after long periods of low precipitation, and due to their slow evolution, it is 
difficult to specify their beginning and end, as well as to estimate their 
magnitude, duration and territorial scope. On the other hand, due to the 
continual growth in water demand for all uses and the threat of climate 
change, in recent years the study of droughts has intensified (Mishra € 
Singh, 2010; Vicente-Serrano, Beguería 8 López-Moreno, 2010). 


Since droughts are a natural phenomenon and threat, their best 
characterization is carried out using indicators that are climatic, hydrological, 
agricultural and/or socio-economic variables, numerically showing the 
potential of droughts to produce deficiencies and damages. A drought index 
is a method based on statistical calculations using indicators. It is aimed at 
quantifying droughts, thereby enabling their comparison among climates and 
dissimilar hydrological regimes. It also helps in the development of correct 
and efficient actions to mitigate the negative impacts of droughts (Mishra € 
Singh, 2010; Fuchs, Svoboda, Wilhite € Hayes, 2014). 


Another important aspect of droughts is their character of multiple time 
scales (multiscaler), which is caused by the delays that occur in the use of 
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water in the different hydraulic systems, so that deficiencies in precipitation 
are first visible in seasonal agriculture, then in  irrigation supplies, 
hydroelectric power generation and other water uses, and finally, in social, 
economic and environmental impacts. According to the above, droughts are 
classified as: meteorological, agricultural, hydrological and socio-economic. 
The latter include environmental damage (Pandey, Sharma, Mishra, Singh, 8 
Agarwal, 2008; Vicente-Serrano et al., 2010; Fuchs et al., 2014). 


Since the end of the last century, drought studies and the development of 
related warning and monitoring systems have been based on the PDSI 
(Palmer Drought Severity Index) and the SPI (Standardized Precipitation 
Index). The former is based on a soil moisture balance (Palmer, 1965; Wells, 
Goddard, € Hayes, 2004) and the latter on a probabilistic approach to 
precipitation (McKee, Doesken, €l Kleist, 1993; Wu, Svoboda, Hayes, Wilhite, 
Ss Wen, 2007). The PDSI has been criticized for not having a multiscalar 
character with respect to time and complicated calculations, besides not 
being sufficiently sensitive to induced changes in precipitation and 
temperature. The SPI has been shown to efficiently use available 
information, while its main criticism is that it uses only monthly precipitation 
data, thereby ignoring other indicators of drought, such as temperature and 
potential evapotranspiration. 


Since droughts are a complex and natural phenomenon with widespread 
impacts, their characterization has recently been sought through indices 
based on multiple variables. Hao and Singh (2015) describe the various 
approaches to constructing multivariate indices, among the best known are 
the use of the soil water balance and latent or hidden variables. A hidden 
variable is constructed with a difference or quotient of variables that have 
great physical significance in droughts, for example, precipitation and 
potential evapotranspiration, and thus the RDI (Reconnaissance Drought 
Index) was proposed with the quotient of these variables (Tsakiris 8 
Vangelis, 2005; Campos-Aranda, 2015), and the SPEI (Standardized 
Precipitation—Evapotranspiration Index) with its difference, which indicates 
conditions of water deficit or excess. The SPEI is a new drought index 
proposed by Vicente-Serrano et al. (2010) and Beguería, Vicente-Serrano, 
Reig and Latorre (2014) whose calculation is similar to the SPI. It is 
multiscalar and sensitive to changes in the historical precipitation and 
temperature records, according to the expectations presented by climate 
change in a zone or region. An application of the SPEI with such an approach 
has been presented by Tórnros and Menzel (2014) for a region in the Middle 
East. 


The objecive of this study is to describe in detail the operating procedure for 
the SPEI and to apply it to the available historical precipitation and average 
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temperature records from the Zacatecas weather station, located in the 
capital city of Zacatecas, in the state of Zacatecas, Mexico. This registry 
covers 86 years in the span of 1930 to 2015. Ten durations of drought 
ranging from 3 to 60 months were processed and three scenarios were 
analyzed for possible climate change. The first considers a progressive and 
linear reduction that reaches 20% of the annual precipitation, the second 
considers a progressive and linear increase that reaches 4*C in the average 
annual temperature, and the third, the most critical, combines the negative 
effects of both types of changes. Based on the percentages defined by SPEI 
for light, moderate, severe and extreme droughts, the results are compared 
and conclusions are formulated. 


Summary of the operational theory 


SPEI calculation 


The SPI was applied using various durations k of droughts, which ranged 
from 3 to 72 months (Vicente-Serrano et al., 2010). The monthly 
precipitation record was processed using moving sums, so that for the 
duration of 3 months, the first sequence, was obtained by adding historical 
data one, two and three, the second sequence was the sum of data two, 
three and four, and so on until the last sequence, which was made up of the 
antepenultimate, penultimate and the last data. Then the number of ns 
sequences that could be formed depended on k and is defined by expression: 


ns = 12:'NA - k+ 1 (1) 


in which NA (for its Spanish initials) is the number of complete years in the 
processed record (> 30 years). The fundamental difference between the 
calculation of SPI and SPEI is that the latter uses as historical data the 
differences (D;i) between monthly precipitation (Pi) and the monthly 
potential evapotranspiration (ETP;;, for its Spanish initials), that is: 
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Ds = Pji= ETPj¡ (2) 


Logically, j varies from 1 to 12 and ¡from 1 to NA. In order to proceed with 
the calculation of the SPEI, first, the sequences of differences designated by 
Df are formed, in which t varies from 1 to ns. When taking into account that 


most of the differences D' are negative, there is a need to use a probabilistic 


model with three fit parameters, whose location parameter u is less than the 
smallest of the sequences to be processed. Vicente-Serrano et al. (2010) and 
Beguería et al. (2014) contrasted four distributions: Log-Normal, Pearson 
type III, Generalized xtreme Value and Log-Logistics; the latter ¡is 
recommend, fitted with the method of biased weighted probability moments 
(Bs), whose equations are (Stedinger, Vogel, € Foufoula-Georgiou, 1993; 
Vicente-Serrano et al., 2010): 


8. = 2) (1-EY-DF with s=0,1,2, (3) 
ns 11 

being: 
10.35 

ES (4) 


In Equation (3), the sequences are in increasing order (DESDILOEDE). 
Beguería et al. (2014) have proposed to use the unbiased Bs when biased 
estimators do not lead to a numerical solution. Stagge, Tallaksen, 
Gudmundsson, Van Loon, € Stahl (2015) have suggested using the General 
Extreme Value distribution, fitted with the maximum-likelihood method. The 
equation for the cumulative probability distribution function [F (x)] of the 
Log-Logistic distribution is (Haktanir, 1991): 


co (5) 


where y > 0,a > O and U < Xmo are the parameters of form, scale and 
location and Xmo is the minimum sequence observed. The values of the fitting 
parameters are estimated with the following expressions (Haktanir, 1991): 
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2-(Bo-3:P2) 
3 6 
(Po-2:P1) ( ) 


ai= (Bo -2-B1) (7) 


y T(1+ y)-T(1-y) 


u=Pp-a T(+y)1M-y) (8) 


in which TF (:) is the gamma factorial function, which was estimated with the 
formula by Stirling (Davis, 1972): 


Msn. 214 O A sl pt) (9) 
12€ 288.82 51840-e€? 2488320-€ 


After calculating the three fit parameters (y, a, u) for each duration k 
analyzed, Equation (5) is applied with x = DF to estimate the non- 
exceedance probabilities F(x) corresponding to each difference. Then, the 
rational numerical approximation is used, developed by Zelen and Severo 
(1972) and presented by Campos-Aranda (2015), to convert F(x) into the 
standardized normal variable Z of zero mean and unit variance, which 
corresponds to SPEI. Light, moderate, severe and extreme droughts are 
defined when the SPEI ranges from O to -1.00, from -1.00 to -1.50, from -— 
1.50 to -2.00 and when it is less than -2.00, respectively. 


Monthly ETP calculation 


Mavromatis (2007) found that the use of simple or complex methods for 
estimating ETP;¡ leads to similar results when applying drought indices such 
as that of Palmer. Based on this result, Vicente-Serrano et al. (2010) 
adopted a simple approach to estimate the ETP);;, through the Thornthwaite 
formula, whose expression is: 
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10-Tt,, 


ETP;; | 


] (10) 


in which Fc is a correction factor function for the latitude of place (LAT) and 
number of days in the month (ndm). Its formula is: 


00 


where N is the maximum sunshine or maximum number of hours with 
average monthly sun. For its estimation in the Mexican Republic, Campos- 
Aranda (2005) developed the following empirical expression: 


N = A + B [sen (30 nm + 83.5)] (12) 


where nm is the number of month, with 1 for January and 12 for December; 
A and B are constant functions of LAT in degrees, with the following 
expressions: 


A = 12.09086 + 0.00266-LAT (13) 
B = 0.2194 - 0.06988:-LAT (14) 


Ttii is the average monthly temperature in *C and IC; is an annual heat 
index, equal to the sum of the 12 monthly indices, which are: 


1.514 
im 22) (15) 
5 
Finally, the exponent m is a function of IC; with the following empirical 


equation: 


m=6.75-107 «IC? -7.71-107 -1C? +1.792-10? -1C,+0.4924 (16) 
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For values of Tt; ¡that are higher than 26.5 *C, there is no influence from 1C;, 
so ETP;¡¡is only a function of Tt;;and is tabulated by Campos-Aranda (2005). 


Hypothetical scenarios of climate change 


Vicente-Serrano et al. (2010) and Ma et al. (2014) have suggested that at 
least two scenarios must be studied as a consequence of climate change 
processes. One is a reduction in precipitation and the other is an increase in 
average temperature. Based on quantitative analyses done by these authors 
and the climatic projections for Mexico by Montero, Martínez, Castillo and 
Espinoza (2010), three scenarios were established for analysis: (1) a 
progressive and linear reduction of 20% in the annual precipitation recorded 
by the historical registry; (2) a progressive and linear increase of 4 *C in the 
average annual recorded temperature and (3) the superimposition of both 
changes in the historical record. The changes mentioned have a direct impact 
on the severity and duration of droughts (Fuchs et a/., 2014). The correction 
of the monthly precipitation record (PM),;) is made based on the following 
equation: 


PM: =PM;(Agots BM se) (17) 


in which Ap is the slope of the reduction and therefore equals the quotient of 
0.20 divided by the number of years NA, and / is the year counter, ranging 
from 1 to NA. The average monthly temperature correction (TM) is carried 
out with equation: 


TM; =TM ;¡ +(Ar -1) (18) 


where Ar is the slope of the increase and therefore equal to the quotient of 4 
“€ divided by NA. 


Processed data and their results 
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Historical records available 


The Zacatecas weather station is located in the city of Zacatecas, which is 
the capital of the state of Zacatecas, Mexico, and which according to 
historical information provided by the Local Office of the National Water 
Commission (Conagua) has operated continuously and has not changed 
location, so the records can be considered reliable. Its geographical 
coordinates are as follows: latitude 22% 45" N, longitude 102% 34' W and 
altitude 2485 m.a.s.!. Its available monthly precipitation (mm) records and 
average temperature (*C), stored in Excel files by CONAGUA of Zacatecas, 
started in January 1953 and are available until December 2015, with missing 
data in April 1986 and several months during the years 2010 to 2013. 


The first missing data were considered equal to the monthly averages and 
the rest was considered equal to the values recorded during the same 
months in the Guadalupe weather station, which is approximately 6 km 
away in a straight line and is located inside the same geographical sub- 
region. This was considered acceptable due to the similarity that both 
records showed at the annual level, both in precipitation and in average 
temperature. On the other hand, in the Climatological Bulletin No. 3 (SARH, 
1980) of the Hydrological Region No. 37 (El Salado), the available records 
from the Zacatecas station for the monthly precipitation and average 
temperature begin in January of 1930 and go up to December of 1978. Then, 
the period from January 1930 to December 1952, without missing data, was 
taken from this Bulletin and with it a joint record of NA = 86 years was 
created. The twelve monthly average values for the integrated precipitation 
record are: 16.2, 9.6, 5.7, 7.5, 17.7, 80.0, 102.1, 97.9, 85.4, 35.4, 12.8 and 
12.0, whose sum is 482.2 mm, which corresponds to the average annual 
precipitation. The average monthly temperatures recorded are: 11.7, 12.7, 
14.8, 17.1, 19.0, 19.0, 17.5, 17.6, 17.0, 16.1, 14.1 and 12.3, with an annual 
average value of 15.7 *C. 


Homogeneity tests of historical records 
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Starting with the monthly historical records, the annual values for 
precipitation and average temperature were integrated, to which three basic 
tests (Helmert, Sequences and Von Neumann) and six specific tests were 
applied: two of persistence (Anderson and Sneyers), two of trend (Kendall 
and Spearman), one of them of change in variability (Bartlett) and another 
of change in the mean, that of Cramer. These tests can be consulted in WMO 
(1971); Campos-Aranda (2005), and Machiwal and Jha (2012). 


The annual precipitation record is perfectly homogeneous, since no general 
or specific test detected any deterministic components. On the other hand, 
the mean annual temperature record was not homogeneous given that the 
basic and specific tests indicated persistence and change in mean. This 
evidence of lack of homogeneity justifies the approach of scenarios related to 
climate change. 


Results from the SPEI with the historical records 


For this study, it was decided to analyze the following ten durations k of 
droughts: 3, 6, 9, 12, 18, 24, 30, 36, 48 and 60 months. In Table 1, the 
results obtained with available historical records are shown for each duration 
according to four concepts: (1) those related to the statistical properties of 
the sequences formed, (2) those corresponding to the fitting of the Log- 
Logistics distribution, (3) those associated with the statistical SPEI indicators 
and (4) those belonging to the drought types that were defined. 


The first group of results provides an idea of the variability and bias of the 
sequences formed with the criterion of moving sums, in such a way that the 
second group, which is the fit parameters of the distribution, is related to this 
behavior. The third group of results is the most important, because it defines 
the quality of the fit obtained and therefore the reliability of the estimates 
made with the SPEI, since its mean, variance and percentage of droughts 
that it defines should be zero, the unit and 50%, respectively. As such 
indicators approximate the values cited, the estimate will be better or more 
reliable. 


As Table 2 shows, the most accurate estimates were those having durations 
of 12, 24 and 36 months, defining approximate percentages of 32.2, 11.5, 
5.0 and 1.3 corresponding to light, moderate, severe and extreme droughts, 
respectively. Figure 1a shows the historical evolution of SPEI with a duration 
of 12 months, observing that wet and dry periods alternate. The longest 
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drought began around sequence 153 (September 1941) and ended in 416 
(August 1963). The date corresponding to a sequence is obtained by clearing 
from Equation (1) the value of NA in years, which is added to the initial year 
of the record. During this period, the most extreme SPEI occurs, with -2.742 
in sequence 321 (October 1954). 
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Figure 1. (a) SPEI evolution of 12 months duration in the Zacatecas weather 
station, state of Zacatecas, Mexico; (b) SPEI evolution (k = 12 months) with 
reduction of 20% in the annual precipitation and increase of 4 *C in the 
average annual temperature and (c) Differences between the SPEI shown (b 
- C). 


SPEI results with modified climate records 


Based on Equations (17) and (18), the historical records were modified to 
incorporate a progressive and linear decrease of 20% in annual precipitation 
and a progressive and linear increase of 4 *C in the average annual 
temperature. Table 2 shows various relevant values for the historical records 
and those for the modified records, according to possible climate change, as 
well as the annual ETP values estimated with the Thornthwaite method. 


Table 1. Characteristic values of the historical climatic records and their 
hypothetical versions according to the probable climate change, in the 
Zacatecas weather station, Mexico. 


historical minimum maximum average 


Annual record of: 
value value 


inicial | final | valor | año | valor | año 
Historical precipitation (mm) 630.9 | 1019.0 | 169.3 | 1969 1019.0 | 2015 482.2 
Average temperature (*C) 15.2 16.2 14.2 | 1985 16.8 1962 15.7 


ETP (mm) with the 736.6 | 763.9 | 704.7 | 1985 | 792.9 | 1962 752.8 
Thornthwaite (THW) method 


Precipitation with progressive 629.4 | 815.2 |145.3|2011 | 880.9 | 1935 432.9 
reduction of 20% (mm) 


Average temperature with 15.2 20.2 15.1 |1931| 20.4 | 2006 17.8 
progressive increase to 4%C 

ETP (mm) of THW and 738.1 | 950.7 | 735.5 | 1931 | 975.3 | 2006 833.6 
average temperature with 

increase 
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Similar to Table 2, in Table 3 the SPEI results from the record modified 
according to the third scenario are presented, /.e., considering the combined 
effect of both changes. Due to space limitation the results of the SPEI with a 
reduction in precipitation and with an increase in average temperature are 
not shown, but they are available from the author. Based on all the SPEI 
results, the duration of 24 months was defined as the ideal for the contrast 
between the values obtained with the historical registry and the modified 
one. This contrast is shown in Table 4. 


Table 2. SPEI results from the monthly historical data from the Zacatecas 
weather station, Mexico. 


N : Duration in months of the analyzed sequences 
umeri 


cal 3 6 9 12 18 24 30 36 48 60 
concept 
s: 


Of 
sequen 
ces 


Number | 1030 | 1027 | 1024 | 1021 1015 1009 1003 997 985 973 
of 
sequen 
ces 


Minimu - = - - - - = - 
m 255.0 389.6 |513.3| 668.6 | 867.1 | 1059. | 1322. | 1487. | 1923. | 2214. 
moving 0 5 3 4 9 
sum 


Maximu - Ml 
m 385.4 | 398.7 | 352.1 | 284.2 | 429.7 | 295.2 | 310.5 43.8 276.6 | 373.2 
moving 
sum 


Arithme | -67.6 - - - - = E Z zS 
tic 135.4 | 203.9 | 273.1 | 412.9 | 552.7 | 692.2 | 832.2 | 1110. | 1389. 
mean 7 4 


Standar | 100.2 
d 135.0 150.4 | 1564 | 214.5 | 234.5 | 278.2 | 289.4 | 337.9 | 380.6 
deviatio 
n 


Asymm | 1.000 
etry 0.799 | 0.640 | 0.548 | 0.525 | 0.539 | 0.546 | 0.451 | 0.234 | 0.145 
coeffici 
ent 


Of the 
FDP 
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(LL3) 


Shape 
parame 
ter (y) 


0.170 


0.148 


0.109 


0.097 


0.080 


0.090 


0.090 


0.084 


0.063 


0.043 


Scale 
parame 
ter (a) 


Locatio 
n 
parame 
ter (u) 


Of SPEI 


Arithme 
tic 
mean 


Varianc 
e 


Minimu 
m 
value 


Number 
of 
negativ 
e 
values 


% of 
drought 
Ss 


52.3 


50.5 


50.6 


48.9 


50.3 


50.8 


50.6 


50.9 


49.4 


Of the 
types 
of 
drought 
Ss 


% of 
light 
drought 
s 


36.2 


32.0 


33.5 


32.8 


29.4 


32.0 


33.6 


33.0 


29.7 


27.9 


% of 
modera 
te 
drought 
Ss 


8.5 


12.2 


10.6 


13.3 


10.7 


17.5 


17.6 


% of 
severe 


6.0 


5.6 


5.8 


4.8 


5.2 


5.7 


4.9 


4.2 


2.1 


3.3 
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drought 
s 


% of 
extrem 


e 
drought 
s 


1.6 


0.8 


1.3 


1.3 


1.0 


1.3 


1.7 


1.9 


1.5 


0.7 


Table 3. SPEI results from the monthly modified records with combined 
effects (reduction in precipitation and increase in average temperature), in 


the Zacatecas weather station, Mexico. 


Numeri 
cal 
concept 
Ss: 


Duration in months of the sequences analyzed 


9 


12 


18 


24 


30 


36 


48 


60 


Of 
sequen 
ces 


Numbe 
rof 
sequen 
ces 


1030 


1027 


1024 


1021 


1015 


1009 


1003 


997 


985 


973 


Minimu 
m 
moving 
sum 


632.2 


766.5 


Maximu 
m 
moving 
sum 


278.9 


165.0. 


395.9 


259.6 


265.5 


Arithm 
etic 
mean 


301.1 


402.3 


Standa 
rd 
deviati 
on 


147.8 


161.3 


225.3 


261.8 


Asymm 
etry 
coeffici 
ent 


0.658 


0.586 


0.671 


0.663 


0.639 


0.596 


0.542 


0.590 


Of the 
FDP 
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(LL3) 


Shape 
parame 
ter (y) 


0.090 


0.076 


0.084 


0.078 


0.074 


0.070 


0.086 


0.104 


Scale 
parame 
ter (a) 


Locatio 
n 
parame 
ter (u) 


Of SPEI 


Arithm 
etic 
mean 


Varianc 
e 


Minimu 
m 
value 


Numbe 
r of 
negativ 
e 
values 


% of 
drought 
Ss 


50.6 


50.2 


48.6 


49.0 


50.1 


49.5 


50.8 


50.2 


49.1 


Of the 
types 
of 
drought 
Ss 


% of 
light 
drought 
s 


34.3 


33.6 


29.7 


30.4 


33.1 


32.3 


33.1 


32.1 


30.8 


% of 
modera 
te 
drought 
Ss 


8.7 


9.6 


12.0 


9.9 


10.9 


10.1 


9.8 


10.2 


12.1 


% of 
severe 


4.9 


5.6 


5.4 


5.3 


5.0 


4.5 


6.2 


7.1 


7.2 


7.3 
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drought 
Ss 


% of 2.7 1.4 1.6 2.0 2.1 1.8 1.5 1.1 1.0 0.3 
extrem 


e 
drought 
s 


Table 4. Number of events and percentages of each drought type obtained 
with an SPEI of 24 months duration from the historical records (Table 2) and 
their contrast with the ones modified according to probable climate change. 


Types of Reduction in Increase in Table 3 
d htand Table 2 
rengnean precipitation | temperature Combined 

minimum value effects 

of SPEI 

Light 323 357 (35.4%) 356 (35.3%) 326 (32.3%) 
(32.0%) > > > 

Moderate 114 90 (8.9%) < 88 (8.7%) < 110 (10.9%) 
(11.3%) < 

Severe 58 (5.7%) 53 (5.3%) < 60 (5.9%) > 45 (4.5%) < 

Extreme 13 (1.3%) 21 (2.1%) > 18 (1.8%) > 18 (1.8%) > 

SPEI -2.334 -2.401 -2.547 -2.453 


The results shown in Table 4 indicate that, in general, light and extreme 
droughts will increase as a result of probable climate change; in contrast, 
moderate and severe droughts will decrease. These results show an anomaly 
in severe droughts corresponding to the record with an increase in 
temperature, which also increases. Regarding the extreme SPEI values, they 
are considered to be in agreement with the modifications applied to the 
historical records. For the third scenario, the most critical, this is shown in 
Table 3 as having more accurate results than those of the 18-month 
duration, with the following percentages for the four types of droughts 
sought: 33.1, 9.9, 5.0 and 2.1. 


Figure 1b shows the evolution of the SPEI for a 12-month duration with the 
record that includes both effects of climate change. It is clear how droughts 
will increase in duration and severity towards the end of the record, 
significantly reducing wet periods. Finally, Figure 1c shows the SPEI 
differences of Figure 1b minus la. 
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Conclusions 


The results from the SPEI based on the historical records from the Zacatecas 
weather station (Table 2), located in the capital of Zacatecas, in the state of 
Zacatecas, Mexico, define approximately 32.2%, 11.5%, 5.0% and 1.3% for 
light, moderate, severe and extreme droughts, respectively. These 
percentages and the following ones were obtained based on drought 
durations of 12 to 30 months. 


In the three climate change scenarios, light and extreme droughts were 
generally found to increase briefly and moderate and severe droughts 
decrease, also slightly in number or percentage. They are defined globally for 
the third scenario (Table 3), the most critical, for combined precipitation 
reduction and average temperature increase, as 32.7%, 10.5%, 4.8% and 
2.0% for each type of drought. 


The evolution graphs of the SPEI, as shown in Figure 1 for a 12-month 
drought duration, allow for accurate visualization and definition of drought 
periods as well as maximum or extreme point values, with respect to both 
the historical record and the records modified according to probable climatic 
change. 


Taking into account the consistency of the results of this study, the 
systematic application of SPEI to the quantification of historical and future 
droughts is recommended, the former based on available historical records 
and the latter on their modified versions, which must consider the possible 
effects of climate change in each zone or region of the country. 


For the climate change scenarios that are considered feasible to occur, SPEI 
results can help to put into perspective the impact on the water balance of a 
region. However, being uncertain how exactly droughts will evolve in the 
future, the SPEI results can only be considered as guidelines. 
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